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Abstract 

The Kuzmin-Toomre family of discs is used to construct potential- 
density pairs that represent flat ring structures in terms of elemen- 
tary functions. Systems composed of two concentric flat rings, a cen- 
tral disc surrounded by one ring and a ring with a centre of attrac- 
tion are also presented. The circular velocity of test particles and 
the epicyclic frequency of small oscillations about circular orbits are 
calculated for these structures. A few examples of three-dimensional 
potential-density pairs of "inflated" flat rings (toroidal mass distribu- 
tions) are presented. 

Keywords: planets: rings - galaxies: kinematics and dynamics. 

1 Introduction 

The presence of flat ring structures is a characteristic feature of the four 
giant planets of the Solar System, in particular, the famous Saturn's rings. 
Planetary rings are the thinnest of all astrophysical discs, with a ratio of the 
thickness to the radius of order 10 -6 pp. On a larger scale, several ring galax- 
ies, or R galaxies, objects in the form of approximate elliptical rings with 
no luminous matter visible in their interiors, are known (see, for instance, 
[21 [3]). Sometimes, as the result of interactions between galaxies, a ring of 
gas and stars is formed and rotates over the poles of a galaxy, originating 
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the so called polar-ring galaxies and references therein). The exact grav- 
itational potential of a ring of zero thickness and constant linear density is 
given by an elliptic integral, which in practice is often approximated by a 
truncated multipolar expansion. 

The potential of a disc with negative mass density whose rim has positive 
mass density can be obtained by a process of complexification of the poten- 
tial of a punctual mass [SI El [7]. This potential was used by Letelier and 
Oliveira [8] to construct superpositions of static bodies with axial symmetry 
in the context of General Relativity. An exact solution representing the su- 
perposition of a disc with a central hole and a Schwarzschild black hole was 
obtained by Lemos and Letelier [9]. The properties of this space-time and 
of space-times that represent superpositions of a Schwarzschild black-hole 
and a ring were investigated in [101 CEH 02] • 

Recently, Letelier |13| constructed several families of flat rings using the 
superposition of Morgan and Morgan [14] discs of different densities. The 
main advantage of these models is that the density and the gravitational 
potential are given in terms of elementary functions. The main purpose of 
this work is to use the same method to construct other potential-density 
pairs representing flat ring structures by superposing members of the classi- 
cal Kuzmin-Toomre family of discs. Because the Kuzmin-Toomre discs have 
infinite extension, so do the resulting rings, while the Letelier rings have a 
finite outer radius or a central hole. All the potential-density pairs can be ex- 
plicitly stated in terms of elementary functions. We construct systems com- 
posed of one and two concentric flat rings and a central Kuzmin-Toomre disc 
surrounded by a flat ring. By using a suitable transformation, some of these 
flat ring systems are used to generate three-dimensional potential-density 
pairs with toroidal mass distributions. Similar toroidal analytical systems 
have been recently obtained by Ciotti and Giampieri [15], and Ciotti and 
Marinacci [16] using a different technique, known as complex-shift method. 

The article is divided as follows. In Section [21 we briefly revise the 
potential-density pair of the Kuzmin-Toomre family of discs. In Section El 
the discs are superposed to generate potential-density pairs that represent 
structures of one and two concentric flat rings, and a central disc surrounded 
by one ring. We calculate the circular velocity of particles moving inside 
these structures and study the stability under radial perturbations of circular 
orbits on the plane of the rings. We also briefly discuss the influence of a 
point mass located at the centre of a ring on the rotation profile and on the 
stability of the ring. In Section 01 we "inflate" some flat ring models by using 
a transformation proposed by Miyamoto and Nagai [T7]. The resulting three- 
dimensional potential-density pairs have toroidal concentrations of matter. 
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In Section [5l we summarize and discuss the results. In Appendices lAltolFl 
we collect, for reference, several equations expressing potential-density pairs, 
circular velocities and epicyclic frequencies for the structures discussed in 
the main text. 

2 Kuzmin-Toomre Discs 

Kuzmin [18] derived a very simple potential-density pair of an axisymmetric 
disc. The potential in cylindrical coordinates (R, z, (p) is given by 

GM 

90 ^/R 2 + (a + \z\y' 

By Poisson's equation, the discontinuous normal derivative on z = intro- 
duces a surface density of mass 



a(R) ' 



2vrG dz 



aM , . 

(2) 



2 =o 27r(i? 2 + a 2 ) 3/2 ' 

The potential (TjQ) is in fact the first member of the Kuzmin-Toomre family 
of discs |19j . The other members can be calculated by using the following 
recurrence relations [201 [2~T] 

<Wi = ^-2^Tl!^ (3) 
ex d 

<?n+l =cr n -o ~~T a -0 "™' n = 0, 1, (4) 

in + 1 aa 

where </>o, o"o is the Kuzmin potential-density pair ([I])-©. The general 
expression for the surface density of the nth-order Kuzmin-Toomre disc is 

(2n + l)a 2n+1 M 

° n ~ 27T( J R2 +a 2 ) n+3/2- ( 5 ) 

For reference the density and potential for the members with n = 1, 2, 3 are 
listed in Appendix lAl 



3 Flat Rings 

3.1 A Family of Single Rings 

In order to generate ring structures we take Kuzmin-Toomre discs with mass 

2-KO?a c . , 

M = -z—rf, 6 

2n+ 1 
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where <r c is a constant with dimensions of surface density. Thus the surface 
density ([5]) is rewritten as 

a c a 2n+3 _ a c 
Gn ~ (R 2 + a 2 )™+ 3 / 2 ~ (l + J R2/a 2 ) n+3/2 ' 

Now let us consider the following superposition 

m 



fc=0 



(1 + i? 2 /a 2 f +m+3/2 ^ 
where = m!/[(m — k)lk\]. Usin. 



m / R 2 \ k 



(9) 



m / R 2 \' 



(10) 



fc=0 

equation (jSJ) takes the form 

= a c (R/a) 2m (n) 
(l+i? 2 /a 2 ) n+m+3/2 ' 

We have that ([lip with m = 1, 2, . . .; n = 0, 1, . . ., defines a family of flat 
rings (characterized by zero density on R = 0) with infinite extension. For 
large R, the density decays as l/i? 2n+3 ; thus in principle, a cut-off radius 
may be defined and this flat ring is considered as finite. The density has a 
maximum at = ( n _^y 2 ) 1//2 . The total mass J\4. < ^ n ' ri ^ of member (m, n) is 
given by 

T r^"' fltlfi d y = r '"', +1)r '" + ) 1 { 2 >^, (i 2 ) 

Jr=o J v =q r(m + n + 3/2) 

where T(x) is the gamma function. As rn, n are integers, ()12p can be further 
simplified to 

(2m + 2n + l)!! v 7 

The potentials associated with (fTTI) can be calculated using a superposi- 
tion with the same coefficients as in ([8]), e.g., the potential associated with 
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m=1,n=0 




Figure 1: The surface density a^ n,n ^ ju c equation (jlip of four members with 
m = 1, 2; n = 0, 1 of the family of flat single rings as function of R/a. 

-( 1 ' ) is 0( 1,0 ) = 0o — 01) etc. The explicit potential-density pairs of four 
members with m = 1, 2; n = 0, 1 are given in Appendix iBl 

In Fig. Q] we show the curves of <j( m ' n ^ ja c as function of R/a for the four 
members of rings listed in equations (134j) - (|37p . The rings with m = 2 have 
larger central holes than those with m = 1. Also note how fast the den- 
sity decays for n = 1 (oc R~ 5 ). These four members have masses .A/f^ 1,0 - 1 = 
4vra c a 2 /3, .A^ 1 ' 1 ) = 4vro- c a 2 /15, 7W (2 ' 0) = 16vra c a 2 /15 andTW^ 2 ' 1 ) = 16vrcr c a 2 /105. 

Some other physical quantities of interest are the circular velocity v c of 
test particles concentric to the flat rings and the epicyclic frequency k of 
small oscillations about the equilibrium circular orbit in the plane of the 
ring. They are calculated with the expressions [22] 

v\ = Rcj> >R , (14) 

« =0,^ + —^-, (15) 

where the subscripts with comma indicate partial derivatives and all expres- 
sions are evaluated on z = 0. The criterion for stability is given by k > 0. 
This study of stability can be considered as an order zero test of stability, 
where the collective behaviour of the particles of the ring is not taken into 
account. 
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Using the potential (|34p . we obtain the following expression for the ro- 
tation curve of the first member of the family of single rings 



v. 



C 1 - ) = R 



2ira c Ga 2 (2R 2 - c 
3(R 2 + a 2 ) 5 /2 



1/2 



(16) 



We have that circular orbits are only possible in the region R/a > V2/2 « 
0.71. Near the hole of the ring, there is not enough mass to produce the 
necessary gravitational force to keep a particle in a circular orbit, however 
this can be compensated if a massive point mass is placed at the centre of 
the ring (see Section 13. 4p . The rotation curves and the expressions for the 
epicyclic frequency «;( m > n ) of the other members are listed in Appendix [Cl 

In Fig. EHa)-(b) we present, respectively, the curves of the circular veloc- 
ity Vc = Vc I V a c Ga and of the epicyclic frequency «( m ' n ) = «;( m > n ) / (a c G/a) x l 2 
as functions of R/a for our four members of rings. We note that the rings 
with m = 2 have larger regions where circular orbits are not possible than 
rings with m = 1 because the former have also larger central holes. The 
possible circular orbits shown are all stable under radial perturbations. In 
Fig. [2](b) we only plotted the curves of epicyclic frequency in the ranges 
where circular orbits are possible. The intervals of stability shown in equa- 
tions (|42 p - (|45p suggest that rings with larger central holes (m = 2) also 
have more unstable orbits. 

3.2 A Family of Double Rings 

Now we consider the following superposition of flat rings, 

(m,n) = cr (m,n)_ 2A ,2 cr (m+l,n-l) A ,4 cr (m+2,n~2) = (R/a) 2 " 1 (l ~ k 2 R 2 /a 2 ) 
(2) (1+ R2/a2) n + m+3/2 ' 

(17) 

where k is a constant. With this superposition, a gap is placed at R/a = 1/k 
on a single ring, thus generating a family of two concentric flat rings. Using 
equation (|13p . the total mass of this superposition is found to be 



. Jmn) -Ka c a 2 2 m+1 m\(2n- 5)\\ r , w . , 9 , 

< ] = ( 2m + 2n + l)!! [(2n " 1)(2n " 3) " 4k (m + 1)(2n " 3) 

+4/c 4 (m + l)(m + 2)]. (18) 

We study in detail the first member of this family with m = 1, n = 2. 
Its non-dimensional surface density is shown in Fig. [3] for the values k 2 = 2, 
k 2 = 1.5 and k 2 = 1. The size of the gap increases for smaller k. 
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(b) 



dm,n) 




Figure 2: (a) The circular velocity vi m ' n ^ = vi" ' /^Ja c Ga equations (f38l) - 
(fSj) and (b) the epicyclic frequency ) / '(a c G I 'a) 1 / 2 equations 

(|42p - (|45p as functions of R/a for the four members with m = 1, 2; n = 0, 1 
of the family of flat single rings. 
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1 2 3 4 5 

R/a 



Figure 3: The surface density <rLj /a c of the first member of a family of 
double rings equation (|17p . as function of R/a, for k 2 = 2, k 2 = 1.5 and 
k 2 = l. 

(I 2) 

The associated potential cp J\ is given by 

= -^^{6x 3 (l+8fc 4 )+x 2 [-16fc 2 (i? 2 +|z| 2 )+a 2 (9-26fe 2 -54^) 

+3a|z|(2 - IQk 2 - 19/c 4 )] - a 2 X [2R 2 {2 + 11A; 2 + 9A: 4 ) + |z| 2 (l - 26A: 2 - 27/c 4 ) 
-a|z|(13 + 82A; 2 + 69k 4 ) - 2a 2 (7 + 28/c 2 + 21/c 4 )] 

-3a 3 (a + \z\){l + A; 2 ) 2 [2# 2 + 7(a + |z|) 2 ] }, (19) 

where x = ^ 2 + (a + M) 2 • Expressions for the circular velocity v^j$ an d 

the epicyclic frequency follow from equation (|19p . and are listed in 

Appendix |Pl 

The curves ^(2) = ^ ( son d nne ) an d K (2) 2 '' = ^ (dashed line) as functions 
of k and R/a are displayed in Fig. SJ Circular orbits are possible in the 
region labeled (A) delimited by the solid curve and the axes; orbits are 
stable in the region (B) delimited by the dashed curve and the axes. Between 
0.58 < k < 1.53 there exist two regions where orbits are unstable and circular 
orbits are not possible: the first begins at the origin and the second is in 
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0.5 1 1.5 

R/a 



(1 2) (1 2) 

Figure 4: The curves v\L = (solid line) and = (dashed line) as 

functions of k and R/a for the first member of a family of double rings. 



an annular region. This is possibly due to the presence of the additional 
gap. We note that most possible circular orbits are also stable, except 
for the interval 0.86 < k < 1.47 where some parts of the double rings 
have unstable circular orbits. Some curves of the dimensionless circular 
velocity = v c(2) I V a cGa and dimensionless epicyclic frequency ftp) = 
K (2) /{ a cG/a) 1 /' 2 for k = 0.5, k = 1 k = 1.5 and k = 2 are depicted in Fig. 
03(a) and Fig. 0(b), respectively. 

3.3 Discs with Flat Rings 

It is also possible to construct potential-density pairs for a Kuzmin-Toomre 
disc surrounded by a concentric flat single ring. For this, we consider the 
superposition 
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(a) 



2.5 



1.5 



v c(2) 



k=0.5 
k=1.0 
k=1.5 
k=2.0 



1 - 




(b) 




Figure 5: (a) The circular velocity v 

,(1,2) 



,(1,2) 



c(2) 



,(1,2) . 
J c(2) ■ 

(b) the epicyclic frequency ^YX' = Km?^ 

first member of a family of double rings as functions of R/a for k 
k = 1 ,k = 1.5 and k = 2. 



/ y 1 o c Ga equation (f46l) and 

k ' /{(JcG/a) 1 / 2 equation (gTJ) for the 

0.5, 
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(2) 

concentrated. The associated potential (jr^ can be expressed as 



Equation (|20p represents the density of a disc with radius R/a = 1/k and a 
flat ring between 1/k < R/a < oo. The total mass of this superposition is 

M?) = 7 w 2 ^ cQ2 w r\(2n-l)(2n-3)-4k 2 (2n-3)+8k 4 ]. (21) 

W (2n+l)(2n-l)(2n-3) LV A ; V ; 

Again we study the first member of this family with n = 2. In Fig. [6] the 
non-dimensional surface density cr|^|/cr c is shown for the values k 2 = 3, 
/c 2 = 2 and k 2 = 1. For large values of fc, the density of the ring is more 
a 

*M = -^g^{x 2 (3-4A; 2 )+ X [a(a+|z|)(3-4fc 2 )+fc 4 (8 J R 2 +8|z| 2 +9ak|)] 
+a 2 [-i? 2 (l + 2fc 2 ) + (a + |z|) 2 (2 + 4A: 2 + 3A: 4 )] }, (22) 
where % = R 2 + ( a + \ z \) 2 ■ By using equation ([22]) . expressions for the 

(2) (2) 

circular velocity irA and the epicyclic frequency can be calculated; 
they are given in Appendix [El 

(2) (2) 

The curves irA = (solid line) and «;W = (dashed line) as functions 
of A; and R/a are shown in Fig. [7J In the region (A), delimited by the solid 
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0.5 



(B) 



1 1 1 1 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

R/a 

— (2) (2) 

Figure 7: The curves v c ^ = (solid lines) and k,J^ = (dashed lines) as 
functions of k and R/a for the first member of a superposition of disc and 
ring. 

curve and the axes, circular orbits are possible and in region (B), delimited 
by the dashed curve and the axes, we have stability. In the limit k = 0, we 
have a pure disc that is stable and all circular orbits are possible. As the 
value of k grows, the influence of the outer ring on the rotation curves and 
stability becomes more pronounced. In the interval < k < 0.86, all circular 
orbits are also stable. When 1.22 < k < 3.21, there exist possible circular 
orbits that are unstable. For k > 3.21, all possible circular orbits are stable. 
Curves for the circular velocity and epicyclic frequency are qualitatively 
similar to those for the two-ring system. 

3.4 A Planet and a Ring 

Now we study a very simple example of a planet located at the centre of a 
ring. Let us consider a point mass with potential eft = —Gm/R and the first 
member of the family of single rings with potential <^( 1,0 ) , equation (|34p . The 
rotation curve and epicyclic frequency of particles orbiting this planet-ring 
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0.4 



0.2 




a=0.04 
a=0.08 
<x=0.12 



2.5 



3.5 



Figure 8: The circular velocity v c = v c /y/a c Ga equation (|23l) for a planet- 
ring system as function of R/a for a = m/(a c a 2 ) = 0.04, a = 0.08 and 
a = 0.12. 



system are given by 



Gm 



Gm 2ira c Ga 2 R 2 (2R 2 - a 



+ 



+ 



R ' 3(R 2 + a 2 f/ 2 
2-KUrGa 2 



1/2 



(2i? 4 + 13i?V - 4a 



R 3 3(i? 2 + a 2 ) 7 /2 



1/2 



(23) 
(24) 



The contribution of the planet to the circular velocity is always positive. As 
was pointed out in Sec. 13.11 circular orbits were not possible near the ring's 
centre due to a lack of mass in this region of space. But if the planet has 
enough mass, the square of the circular velocity can be made non-negative 
everywhere. This happens if m/(a c a 2 ) > 47T\/5/375 ~ 0.075. In Fig. [8] we 
display some rotation curves, where v c = v c j \J a c Ga is the dimensionless 
velocity and where we defined a = m/(a c a 2 ). By equation (|24p . the contri- 
bution of the central mass to the epicyclic frequency is positive, so it helps 
to stabilize the circular orbits in the plane of the ring. In fact, it can be 
shown that if 



m 



a r a* 



> ^r(65Vl9 - 41^)^/70 - 2^/665 « 0.15, 



3675 



(25) 
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2 

R/a 



Figure 9: Isodensity curves of the mass density p^ 1 ' ' / (a c / a) equation (|27l) 
for a thick single ring as functions of R/a and z/a. Parameters: (a) b/a = 1 
and (b) b/a = 0.5. 



the epicyclic frequency is always non-negative. 



4 Thick Rings 

Until now all ring structures presented were flat, i.e., had infinitesimal 
thickness. One way to generate three-dimensional potential-density pairs 
from these flat ring models is to employ the same transformation used by 
Miyamoto and Nagai [T7] to "inflate" the Kuzmin-Toomre discs. 

As a first example, we take the potential c^ 1,0 ) of the first member of the 
family of single rings and apply a transformation \z\ — ► V z 2 + b 2 , where b 
is a non-negative constant, to equation (|34p . The resulting mass density is 



calculated directly from the Poisson equation in cylindrical coordinates 



1 



4vrG 



\RR + 



R 



+ 



(26) 



We obtain 

(i,o) 



a c a 2 b 2 



P 



2?[R 2 + (a + 2 ] V2 



{ai? 4 +i? 2 (a+0 [a 2 + 2£(3a + £)] +2£ 2 (a+0 3 }, 

(27) 

where £ = V z 2 + b 2 . 

By inspection, the mass density equation (|27p is non-negative and free of 
singularities. Fig. [I2a)-(b) show some isodensity curves of p^ 1 ' ' / (a c /a) as 
functions of R/a and z/a with parameters b/a = 1 in Fig.[9](a) and b/a = 0.5 
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in Fig.[9jb). For this last parameter value, it is seen that we have a toroidal 
mass distribution. We also "inflated" the single ring with potential c^ 1,1 ) 
and found a three-dimensional mass density with similar properties. 

It is also interesting to compute the rotation curve and epicyclic fre- 
quency of particles orbiting the symmetry plane z = of the "inflated" 
ring. We get 



J c(th.) 



R< 



2-ko c Gcl 2 [2R 2 - (a + 6) (a - 2b)] 



1/2 



5/2 



,(1,0) 



2-ko c Gcl 2 



R 2 + (a + bf 

^ [2R A + R 2 (a + b) (13a + 106) 



(28) 



3 [R 2 + (a + 6) 

-4(a + 6) 3 (a-26)]} 1/2 



(29) 



where the subscript th. refers to thick structures. As the thickening parame- 
ter 6/a is increased, the regions with impossible and unstable circular orbits 
shrink until they disappear when 6/a > 1/2. On the other hand, Fig. sug- 
gest that prominent toroidal mass distributions are formed for 6/a < 1/2, 
and these still have regions with impossible and unstable circular orbits. 

If we apply a Miyamoto-Nagai transformation to the potential of two 
concentric flat rings, we expect to generate two toroidal mass distributions. 
For instance, the potential 0(2) equation (fT9j) yields a somewhat involved 
expression for the density, 



.(1,2) 
P{2) 



a c a 2 b 2 



■ ^5ak 4 R 8 + R 6 [2f(3 - 8k 2 + 24k 4 ) 



70e[R 2 + (a + O 2 ] U/2 
+35a£; 4 £(lla + lOf ) + 35£; 2 a 3 (£; 2 - 2)] + i? 4 (a + f ) [6£ 4 (3 - 8k 2 + 24k 4 ) 

+a^ 3 (48-128A; 2 +699A: 4 )+5a 2 ^ 2 (7-112A; 2 +126/c 4 )-35a 3 ^(l+20A; 2 )+35a 4 (l-2A; 2 

+i? 2 (a+e) 3 [6e 4 (3-8/c 2 +24A: 4 )+2a^ 3 (39-104A; 2 + 172A; 4 )+5a 2 C 2 (21-112A; 2 ) 

+5a 3 (7a+56£)] +2£ 2 (a+£) 5 [£ 2 (3-8£; 2 +24£; 4 )+2a£(9-24£; 2 +2£; 4 )+35a 2 ] J, 

(30) 



where £ = \J z 2 + b 2 . 

(1 2) 

Some isodensity curves of the mass density p\ 2 l /(a c /a) are displayed in 
Fig. HoTa)-(b) with parameters k = 1, b/a = 0.5 in Fig. HOT a) and k = 1, 
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Figure 10: Isodensity curves of the mass density pLx /(a c /a) equation ([30]) 
for a thick double ring as functions of R/a and z/a. Parameters (a) k = 1, 
b/a = 0.5 and (b) k = 1, b/a = 0.3. 



b/a = 0.3 in Fig. fTOTb). For this set of parameters, the mass density is 
positive everywhere. We also tested it for some other values of k and b/a 
and did not find regions with negative mass density or singularities. 

For this two-toroid system, we also computed the rotation curve and 
epicyclic frequency. The explicit expressions are given in Appendix [F] Fig. 
frfT a) shows some zero velocity curves of equation ([50]) and Fig. [TU b) curves 
of zero epicyclic frequency equation (|51 [) as functions of R/a and k for the 
thickening parameters b/a = 0, b/a = 0.25 and b/a = 0.5. Here, thickening 
also helps to stabilize circular orbits and shrinks the regions with negative 
square velocities. 



5 Discussion 



We presented some families of analytical potential-density pairs that repre- 
sent one or two concentric flat rings that were constructed by superposing 
members of the Kuzmin-Toomre family of discs, and also a disc surrounded 
by a flat ring. All potential-density pairs can be expressed in terms of ele- 
mentary functions. Although the models have infinite extension, the surface 
density decays with the distance very fast, and, in principle, one can insert 
a clear cut-off and consider them as finite. For these models, we studied 
the circular orbits of test particles concentric to the rings and made a first 
study of the stability of these structures by radial perturbations of circular 
orbits (epicyclic frequency). In general the two-ring structures tend to be 
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more unstable than the single rings, and a surrounding ring also helps to 
desestabilize a disc. On the other hand, if we put a planet located at the 
centre of a ring it tends to increase stability in the radial direction. 

We also presented two examples of "inflated" rings that were obtained 
by applying a Miyamoto-Nagai like transformation on the above-mentioned 
single and double flat rings. The resulting three-dimensional mass densi- 
ties have a toroidal shape, are non-negative and free of singularities. The 
thickening of flat ring structures increases stability in the radial direction. 

As was mentioned, a first stability test of the presented ring structures 
was based on the study of stability of individual particles. The study of 
stability of orbits perpendicular to the plane of flat rings is not trivial because 
of the abrupt change of matter density as a particle crosses the plane of the 
ring. For instance, Saa and Venegeroles [23] showed that the presence of a 
thin plane destroys the integrability of orbits of particles that cross the plane 
and these orbits tend to be chaotic. This effect is similar to the behaviour 
of the Chua circuit, a well-known case of chaotic system [23]. Recently, [25] 
studied the stability of orbits crossing Morgan and Morgan discs, thin finite 
dies. 

A more realistic approach should take into account collective effects. The 
stability of the gravitating system will be determined by the behaviour of 
the perturbation of the matter content of the disk. Some examples of a 
General Relativistic stability analysis of discs can be found in [261 [27] . 
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A Potential-density pairs for Kuzmin-Toomre discs 

The explicit expressions for the density and potentials for Kuzmin-Toomre 
discs with with n= 1, 2, 3 are 

(31) 

5a 5 M , GM r 9 . . ,. 91 

2vr(i? 2 + a 2 ) 7 / 2 3 [i? 2 + (a + |z|) 2 ] 

x [i? 2 + (a + |z|) 2 + o(a + |z|)] - a 2 [R 2 - 2(a + |z|) 2 ] }, (32) 
7a 7 M GM 



03 



^ = ~ r o , T^FH + (° + I 2 !)'] 



2vr( J R 2 + a 2 )9/ 2 ' 5 [R 2 + (a + N)2 ] V 2 

+a[,R 2 + (a + |z|) 2 ] 2 (4a + 5|z|) + a 2 [,R 2 + (a + |z|) 2 ] [-i? 2 + 5|z|(a + \z\)] 
+a 3 (a + |z|)[2 J R 2 + 7(a + |z|) 2 ]}. (33) 
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B Potential-density pairs for flat one rings 



The explicit expressions for the density and potentials of four members of 
flat one rings with m = 1, 2; n = 0, 1 read 



a 



(1,0) 



a c {R/aY 



r (l,l) 



(l + fl 2 /a 2 ) 5/2 ' 
+a 2 + 3a\z\) , 

<Jc {R/af 
" (l + i? 2 /a 2 ) 7/2 ' 



r (2,0) 



a c (R/a) 



a 



(2,1) 



(l + i? 2 /a 2 ) 7/2 ' 
8i? 2 + 8|z| 2 + 9a| 
^ c (i?/a) 4 



\2/ D 2 



0( 2 '°) = - 



2-Ko c Ga? . 


3[i?2 + ( a+ | z |)2 


3/2 V 


27rcr c Ga 2 




15[i? 2 + (a + |z|) 


2] 5 /2 L 


z| 2 + a 2 + 3a|z|) 




27rcr c Ga 2 




15[i? 2 + (a+|z|) 


2j 5 / 2 i 


M) 2 }, 




27T(T c Ga 2 





" (l + i? 2 /a 2 ) 9/2 ' " 105[ J R 2 + (a + |z|) 2 ] 7/2 

x(8R 2 + 8|z| 2 + 13a 2 + 24a\z\) + a 2 [R 2 + (a + |z|) 2 ] 

x(lli? 2 - 13|z| 2 -28a 2 -41a\z\) + 3a 3 (a+ |z|) [2i? 2 + 7(a + \z\) 2 



(34) 

i? 2 (2i? 2 + 2|z| 2 
(35) 

[i? 2 + (a + |z|) 2 ] 
(36) 

{[R 2 + (a + \z\) 2 ] 2 



(37) 
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C Rotation curves and epicyclic frequency for flat 
one rings 

The expressions for the circular velocity v^' n ^ and the regions with non- 



negative velocity for the members of flat one rings (|34p ~ (|37p are given by 



,(1,0) 



R 



2ira c Ga 2 (2R 2 - a 2 ) 
3(R 2 + a 2 ) 5 / 2 

2ira r Ga 2 



1/2 



„(1,0) 



R a/2 
> for - > — ss 0.71, (38) 
a 2 



1 5(^ + ^ (M4 + 13BV - 4 ° 4: 

2 



1/2 



;P<°) = R 



a 

2ira c Ga 2 
15(R 2 + a 2 ) 7 / 2 



0.54, 



(39) 



(8i? 4 - 8R 2 a 2 - a 4 ) 



1/2 



,P.°> > for * > ^[+2 « 1.05, 
a 2 



^ = R 



2ira r Ga 



W5(R 2 + a 2 ) 9 / 2 



(8i? b + 72i?V - 45#V - 4o b ) 



i? 



1/2 



A 2 ' 1 ) > for - > 0.81 



(40) 



(41) 
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The corresponding expressions for the epicyclic frequency K^ m ' n ^ and the 
regions of stability read 



-(i,o) 



2ira c Ga 2 
3(R 2 + a 2 ) 7 / 2 



(2# 4 + 13i?V - 4a 



1/2 



stable for * > vV^ZH w 
a 2 



-(1,1) 



27T<7 c Ga 2 ,, ir)(i , ,, r ,i 2 



15(fl 2 + a 2 ) 9 / 2 



(2fT + 3i? 4 a^ + 90iTa 4 - 16a 
i? 



1/2 



r (2,0) 



stable for — > 0.42, 
a 

2ira r Ga 2 4 2 iro- 1 



15(i? 2 + a 2 ) 9 / 2 



(8i2° + 72irV - 45fTa 4 - 4a D 
R 



1/2 



.-(2,1) 



stable for — > 0.81, 
a 

2ira c Ga 2 fn n8 , n 6 2 ( ^ 11d4 ^s 



(8i? s + 8i?V + 711# 4 a 4 - 250i? 2 a 6 - 16a 



105(i? 2 + a 2 ) n /2 



i? 



stable for — > 0.64. 



(42) 
(43) 
(44) 

nl/2 

5 

(45) 



D Rotation curves and epicyclic frequency for flat 
two rings 

(1 2) 

The expressions for the circular velocity v\L an d the epicyclic frequency 

K}„\ for a particular member of two flat rings studied in Section 13.21 are 
given by 



v {1 > 2) - R 

c(2) ~ JX 



2ira r Ga 2 



[2R 6 (3 - 8k 2 + 24k 4 ) - 3i? 4 a 2 (-ll + 48A; 2 + 24fc 4 ) 

(46) 



W5(R 2 + a 2 ) 9 / 2 

+ 18i? 2 a 4 (6 + 5k 2 - k A ) - a 6 (24 - 8A; 2 + 3k 4 )] } 1/2 



k {1 ' 2) - { 



2i\a c Ga 2 



■ [2i2 8 (3 - 8A: 2 + 24k 4 ) + i? 6 a 2 (27 - 16fc 2 + 552& 4 



ll05(i? 2 + a 2 )"/ 2 

6i? 4 a 4 (10 + 237A; 2 + 87A; 4 ) + J R 2 a 6 (768 + 500A; 2 - 93/c 4 ) 

-4a 8 (24-8fc 2 + 3fc 4 )]} 1/2 . 



(47) 
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E Rotation curves and epicyclic frequency for discs 
with flat rings 

(2) 

The expressions for the circular velocity trA and the epicyclic frequency 

(2) 

k)J, for a particular member of the disc with a fiat ring studied in Section 
13.31 are given by 

v ^ = R { lEwf^W [RHs " 4 " 2 + 8 " 4) " 2i?V( ~ 6 + 13k2 + 4 " 4) 

+a 4 (24 + 8fc 2 -k A )]} 1 ' 2 , (48) 

K S = { 15(^)9/2 ^ - ^ + + - ^ + 12 ^) 

-45i? 2 a 4 A: 2 (4 + k 2 ) + 4a 6 (24 + 8/fc 2 - A; 4 )] } 1/2 . (49) 

F Rotation curves and epicyclic frequency for thick 
two rings 

The expressions for the circular velocity w^Hh anc ^ ^ e e Pi cvcuc frequency 
K (2)th ^ or a particular member of thick two rings studied in Section U] are 
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given by 

»gh = A 2nacGa2 9/2 {2^ 6 (3 - 8fc 2 + 24^) 

c(2)th - I 105 [R 2 + (a + 6) 2 1 9/2 1 



+3# 4 [-a 2 (-ll + 48/c 2 + 24/c 4 ) + 3a5(6 - 16/c 2 + 13/c 4 ) + 26 2 (3 - 8/c 2 + 24/c 4 )] 
+3i? 2 (a + b) [6a 3 (6 + 5A; 2 - k 4 ) + 26 3 (3 - 8A; 2 + 24A; 4 ) + 10a& 2 (3 - 8£i 2 + 3k 4 ) 
-a 2 6(-45 + 64A; 2 + 39fc 4 )] - (a + bf [a 3 (24 - 8A; 2 + 3A; 4 ) - 26 3 (3 - 8A: 2 + 24fc 4 ) 

+3a6 2 (-12 + 32A: 2 + 9k 4 ) - 6a 2 fe(ll + 8k 2 - 3k 4 )] }} 1/2 , (50) 

j 2Tra r Ga 2 



-{2i? 8 (3-8fc 2 + 24/c 4 



W \l05[i?2 + (a + b) 2 ] 11 ' 2 
+R 6 [a 2 (27 - Wk 2 + 552/c 4 ) + a&(66 - 176/c 2 + 843A; 4 ) + 146 2 (3 - 8k 2 + 24/c 4 )] 
+3i? 4 (a + 6) [-2a 3 (10 + 237A: 2 + 87A: 4 ) + a 2 6(97 - 576A; 2 + 237A; 4 ) 
+2a6 2 (51 - 136/c 2 + 303/c 4 ) + 10o 3 (3 - 8k 2 + 24k 4 )] + R 2 (a + b) 3 
x [a 3 (768 + 500/c 2 - 93/c 4 ) - 12a 2 6(-40 + 116/c 2 + 51k 4 ) 
+15a6 2 (24 - 64/c 2 + 45/c 4 ) + 266 3 (3 - 8/c 2 + 24/c 4 )] - 4(a + 6) 5 
x [a 3 (24 - 8/c 2 + 3k 4 ) + 6a 2 6(-ll - 8/c 2 + 3k 4 ) + 3a6 2 (-12 + 32/c 2 + 9k 4 ) 

-26 3 (3 - 8A: 2 + 24fc 4 )]}} 1/2 . (51) 
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